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1.0 


INTRODUCTION 


This  annual  report  for  the  contract  No.  F49620-83-C-0060  to  per¬ 
form  droplet  measurement  research  covers  the  period  15  January  1983  to 
15  January  1984. 

During  this  first  year  of  the  contract  considerable  progress  was 
experienced  as  indicated  in  the  Statement  of  Work  (Section  2.0).  Two 
techniques  referred  to  as  V/I  and  IMAX  were  studied  and  since  the  two 
are  reasonably  independent  the  technical  sections  are  divided  to  clearly 
identify  the  progress,  achievements  and  limitations  associated  with  each 
technique. 
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2.0 


STATEMENT  OF  WORK 


The  statement  of  work  as  shown  in  the  original  proposal  and  con¬ 
tract  is  included  here  for  the  three  years  of  the  contract.  Next  to 
each  task  there  is  an  asterisk  (*)  or  a  triangle  (A),  the  asterisk  indi¬ 
cating  completion  of  the  task  and  the  triangle  indicating  that  some  work 
was  done  during  this  year. 


Year  1  -  Droplet  Measurement  Research  (Basic  Research  Element) 


*Task  1.0  -  Evaluation  and  assessment  of  previous  results  to  provide  a 
basis  of  departure  for  the  present  research. 

*Task  2.0  -  Plan  the  research  to  study  two  advanced  droplet  sizing  con¬ 
cepts  (V/I  and  IMAX)  both  analytically  and  experimentally. 

*Task  3.0  -  Theoretical  definition  of  the  two  concepts  based  on  funda¬ 
mental  scattering  theories. 

*Task  4.0  -  Breadboard  experimental  setups  to  test  and  study  both  tech¬ 
niques  . 

*Task  5.0  -  Experimentally  determine  the  intensity  of  the  light  scat¬ 
tered  by  droplets  immersed  in  a  Gaussian  beam.  Compare 
these  results  with  present  scattering  theories  and  correla¬ 
tions. 

*Task  6.0  -  Explore  signal  characteristics  of  the  two  techniques  under 
ideal  conditions,  and  establish  optimum  regions  of  valid¬ 
ity. 

ATask  7.0  -  Analyze  above  data  and  compare  to  the  theoretical  descrip¬ 
tion  obtained  in  Task  3.0,  to  determine  optimum  selection 
of  parameters. 

*Task  8.0  -  Report  the  results  of  each  of  the  program  tasks  and  develop 
a  publication  for  submittal  to  a  referee  journal. 
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2.2 


Task  1.0 

Task  2.0 
ATask  3.0 
*Task  4.0 
Task  5.0 
Task  6.0 
*Task  7.0 
Task  8.0 

ATask  9.0 


Evaluate  and  assess  previous  year's  results  and  recent 
relevant  developments  and  re-establish  the  direction  of  the 
proposed  investigation. 

Plan  the  research  for  the  second  year  based  on  information 
obtained  to  this  point. 

Formulate  an  analytical  model  to  compute  the  optical  probe 
volume  and  check  experimentally  verify  the  result. 

Explore  the  signal  characteristics  of  a  monodispersed  spray 
when  laser  beams  are  interfered  by  a  real  spray. 

Evaluate  how  regions  of  established  validity  under  ideal 
conditions  are  affected  by  beam  blockage. 

Formulate  statistical  analysis  to  predict  the  effect  of 
spray  blockage  on  signal  characteristics. 

Define  and  establish  a  processing  method  for  the  concepts 
studied.  Perform  preliminary  spray  measurements. 

Analyze  above  data  and  compare  to  the  statistical  model 
obtained  in  Task  6.0  to  reassess  the  validity  of  the  opti¬ 
mum  parameters  obtained  under  ideal  conditions. 

Report  the  results  of  each  of  the  program  tasks  and  develop 
a  publication  for  submittal  to  a  referee  journal. 


2.3  Year  3  -  Droplet  Measurement  Research  (Basic  Research  Element) 


Task  1.0  - 

Task  2.0  - 
ATask  3.0  - 

ATask  4.0  - 

ATask  5.0  - 


Evaluate  and  assess  developments  and  the  results  of  the 
first  two  years  and  redefine  the  critical  aspects  of  the 
research. 

Plan  the  research  for  the  third  year  based  on  available 
information.  Consider  alternative  methods  or  new  concepts. 

Perform  measurements  in  real  sprays  and  explore  the  signal 
characteristics  and  establish  the  limitations  of  the  tech¬ 
niques  . 

Test  processing  methods  chosen  during  second  year,  in  the 
presence  of  real  sprays,  and  develop  corrections  if  neces¬ 
sary. 

Compare  the  results  obtained  with  the  various  methods  for 
the  same  spray  under  various  spray  conditions. 
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Task  6.0 
Task  7.0 


Analyze  above  data  in  conjunction  with  statistical  models 
to  establish  the  limitations  of  the  techniques. 

Report  the  results  of  the  entire  program  and  develop  a  pub¬ 
lication  for  submittal  to  a  referee  journal. 
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3.0 


ASSOCIATED  RESEARCH 


3 . I  Research  Associated  with  the  Visibility/Intensity  Technique 

The  interferometric  technique  to  measure  the  size  and  velocity  of 
particles  flowing  in  a  fluid  has  been  in  existence  for  approximately  one 
decade.  Farmer*  introduced  the  concept  of  using  visibility  for  particle 
sizing  in  the  forward  scattering  direction.  Since  then,  many 

researchers  have  contributed  to  the  development  of  the  technique. 

2 

Robinson  and  Chu  ,  for  instance,  rederived  Farmer's  results  using  a  more 
rigorous  approach,  and  the  definition  of  visibility  through  the  first 
order  Bessel  function  was  confirmed  experimentally  in  forward  scatter. 
However,  contrary  to  what  was  through  at  the  time,  this  visibility  rela¬ 
tionship  could  not  be  used  at  any  angle  other  than  the  forward  or  zero 
degree.  Adrian  and  Orloff  showed  that  it  could  not  be  used  in  the 
backscatter  mode  as  predicted  by  Mie  scattering  theory.  This  finding 
was  also  confirmed  by  Roberds^,  who  in  addition  showed  that  the  config¬ 
uration  of  the  receiving  optics  would  also  affect  the  relationship 
between  visibility  and  size. 

The  biggest  known  limitation  of  the  visibility  technique  was  its 
inability  to  measure  particles  in  a  dense  field.  The  main  reason  being 
that  the  probe  volume  when  observed  in  the  forward  direction  is  big  and, 
therefore,  the  criterion  of  single  particle  measurement  is  not  always 
met.  In  addition,  the  technique  dictated  that  the  largest  measurable 
droplet  was  of  the  order  of  the  fringe  spacing.  This  imposes  a  limita¬ 
tion  of  about  200  pm  to  the  largest  measurable  particle  and  also  results 


in  a  large  probe  volume,  since  a  minimum  number  of  fringes  is  needed 
(typically  8)  to  accurately  process  a  Doppler  signal. 

Bachalo^  formulated  a  mathematical  model  based  on  classical 
optics  where  the  Mie  scattering  is  approximated  with  refraction,  dif¬ 
fraction  and  reflection.  That  model  was  used  to  predict  the  visibility 
of  spherical  particles  in  off-axis  directions.  Experimental  results 
using  raonodispersed  droplets  confirmed  the  validity  of  the  model  under 
many  conditions.  Pendleton^  later  confirmed  these  results  using  a 
sophisticated  numerical  model  that  utilizes  the  full  Mie  equations. 
Recent  work  conducted  under  this  contract  at  Spectron  Development  Labor¬ 
atories  shows  that  there  are  two  major  limitations  on  the  existing 
models.  First,  all  the  theories  developed  thus  far  predict  the  Mie 
scattering  from  beams  of  uniform  intensity.  Second,  the  formation  of 
the  fringe  pattern  at  the  probe  volume  may  be  affected  by  secondary 
scattering  of  particles  immersed  in  the  laser  beams  before  they  cross. 

As  a  result  of  these  limitations,  the  relationship  between  visibility 
and  size  is  not  straight  forward  and  considerable  errors  can  be  made.  A 
model  proposed  by  M.  L.  Yeoman  et  al^  pursues  the  solutions  of  the  above 
problems.  It  consists  of  two  concentric  interferometric  patterns  of  two 
colors  such  that  the  small  one  establishes  a  region  of  the  big  one  where 
the  hyperbolic  variation  of  the  visibility  is  negligible.  Also  digital 
analysis  of  the  signals  helps  establish  some  that  are  in  error.  The 
results  reported  here  show  that  using  the  intensity  of  the  pedestal  of 
the  scattered  light  in  addition  to  the  visibility  will  eliminate  many  of 
these  errors,  significantly  improving  measurement  accuracy. 
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3.2  Research  Associated  with  the  I MAX  Technique 

Several  nonintrusive  single  particle  counters  using  absolute 
scattered  light  have  been  proposed  and  they  can  be  divided  into  two 
groups.  In  the  first  group^-*^  Gaussian  laser  beams  are  used  and  var¬ 
ious  mathematical  inversion  techniques  have  been  proposed  to  extract  the 
true  size  distribution  from  the  signal  amplitude  distribution.  Several 
assumptions  are  needed  to  perform  such  inversion,  such  as  assuming  the 
form  of  the  size  distribution^^  or  a  uniform  average  velocity  for  all 
the  particles  regardless  of  size^^.  A  method  to  account  for  pinhole 
masking  has  also  been  described^ The  second  group  of  NSPC's  uses  a 
probe  volume  of  known  illumination  to  avoid  the  ambiguity  of  particle 
position.  The  technique  described  here  falls  in  this  group.  Several 
approaches  have  been  proposed^* ^>13,14)  to  accomplish  this  task. 
Apertures  are  used  to  cut  off  the  edges  of  the  Gaussian  beam  and  thus 

Image  a  probe  volume  of  almost  uniform  illumination''  *  .  One  of  the 

r  >  -!  \ 

apertures'  had,  In  addition,  two  wedges  that  when  properly  oriented 
can  be  used  to  obtain  the  velocity  vector  in  two  components  (magnitude 
and  direction).  Ring-shaped  probe  volumes  have  also  been  produced^1^ 
by  a  TEMqi  laser  and  a  ring-shaped  aperture.  The  measurement  volume  is 
defined  by  the  intersection  of  apertures  in  front  of  two  photomulti¬ 
pliers  . 


-7- 


84-2286-06/36 


4.0 


THE  VISIBILITY /INTENSITY  TECHNIQUE 


In  this  section  we  describe  some  of  the  problems  associated  with 
the  visibility  technique,  and  how  many  of  these  problems  can  be  cor¬ 
rected  or  avoided  combining  the  visibility  with  the  intensity  of  the 
scattered  light. 

4.1  Errors  in  the  Visibility  Measurements 

The  prediction  of  the  droplet  size  based  on  the  visibility  of  the 
light  scattered  by  the  droplet  crossing  an  interferometric  pattern  of 
fringes  is  dependent  on  the  pattern  itself.  For  a  perfect  system  the 
fringe  visibility  should  be  100%  and  the  light  scattered  should  be  pri¬ 
marily  of  one  type.  There  are,  however,  many  reasons  in  practical 
environments  that  contribute  to  alter  the  fringe  pattern  and,  therefore, 
confuse  the  relationship  between  visibility  and  size.  The  result  is 
that  apparent  sizes  instead  of  true  sizes  are  often  measured.  As  an 
example,  droplets  interacting  with  the  laser  beams  before  they  cross 
will  randomly  reduce  (or  even  destroy)  the  fringe  visibility.  If  a 
droplet  is  measured  at  such  time,  it  will  appear  bigger  than  its  true 
size.  There  are  several  factors  that  produce  an  error  in  the  visibility 
measurement  and  they  are  listed  below.  They  are  divided  into  two 
groups:  the  ones  that  reduce  the  visibility  and  the  ones  that  increase 
it. 

A)  Factors  that  Reduce  the  Visibility 

1)  Particles  prior  to  crossover  destroy  fringe  con¬ 
trast.  This  is  a  function  of  the  spray  density. 
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2)  Beam  excursl  is  due  to  turbulent  hot  media. 

3)  Hyperboloidal  reduction  of  visibility. 

4)  Reflection/refraction  can  add  destructively  for 
large  droplets  immersed  in  a  nonuniform  (typically 
Gaussian)  beam. 

5)  Multiple  particles  in  probe  volume. 

6)  Increase  in  signal  background.  Normally  from  many 
little  drops  present  in  the  probe  volume. 

7)  Loss  of  fringe  contrast  due  to  faulty  components 
(beamsplitter,  etc.) 

B)  Factors  that  Increase  the  Visibility 

1)  Out  of  focus  drops  "masked”  by  the  pinhole. 

2)  Low  level  signals  in  the  presence  of  a  noisy  back¬ 
ground. 

3)  High  time  rise  on  leading  edge  of  signal. 


4 .2  Description  of  the  V/I  Technique 

This  technique  makes  use  not  only  of  the  visibility  of  the  signal 
but  also  the  peak  intensity  of  the  pedestal.  Both  parameters  are  avail¬ 
able  in  the  signal  and  their  cross-correlation  can  be  used  to  eliminate 
faulty  signals  produced  in  many  practical  environments.  This  technique 
will  especially  prevent  small  particles  (high  visibility)  from  appearing 
as  large.  Figure  4.1  shows  a  Doppler  trace  where  both  the  visibility 
and  the  peak  intensity  of  the  pedestal  are  defined.  There  is  a  correla¬ 
tion  between  the  size  of  the  droplet  and  the  amount  of  scattered  light 
given  by  Mie  theory.  This  correlation  can  be  used  to  eliminate  signals 
with  an  apparently  different  size.  The  proposed  logic  is  as  follows: 
droplets  that  produce  a  certain  visibility  must  have  a  given  size; 
hence,  they  must  scatter  light  with  a  given  intensity  (characterized  by 
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I 


Figure 


Ip)*  Two  exceptions  are  contemplated:  first,’  droplets  with  the  correct 
visibility  will  scatter  different  amounts  of  light  due  to  the  Gaussian 
nature  of  the  probe  volume's  intensity;  second,  droplets  with  an  erro¬ 
neous  visibility  will  not  scatter  light  with  an  intensity  corresponding 
to  their  apparent  size. 

The  V/I  method  will  then  establish  intensity  limits  for  every 
measured  visibility.  This  will  produce  a  well  established  probe  volume 
as  a  function  of  size  and  will  reject  droplets  measured  with  an  erro¬ 
neous  visibility.  For  instance,  assume  that  the  size  corresponding  to  a 


visibility  is  dp  Then  the  light  scattered  will  have  a  pedestal  with 
maximum  intensity  Ip^(max).  Since  this  intensity  would  only  be  detected 
from  droplets  crossing  through  the  middle  of  the  probe  volume,  more 


relaxing  limits  are  proposed.  That  is  Ip  will  be  accepted  when 
Xp^(min)  Ip^(min)  <  Ip^  <  Ip^(max).  Obviously,  the  broader  the  limits, 
the  larger  the  error  that  can  be  allowed  with  the  visibility  technique. 


Figure  4.2  shows  three  cases  that  illustrate  the  correcting  prop¬ 


erties  of  the  V/I  method.  Figure  4.2a  shows  a  droplet  of  size  with 
visibility  and  intensity  I—*-- .Since  this  last  one  is  within  estab¬ 
lished  limits,  the  signal  is  accepted  as  valid.  Figure  4,2b  shows  the 


signal  produced  by  a  droplet  of  diameter  but  the  visibility  was 
reduced  to  due  to  fringe  contrast  reduction.  Since  the  visibility  is 
Vp  this  signal  would  be  interpreted  as  having  a  diameter  dp  However, 


the  peak  intensity  of  the  pedestal  lp  corresponds  to  its  true  size  (d2 ) 
and  is  smaller  than  Ip^(min)  (to  which  it  is  compared)  and  hence,  the 


signal  is  rejected.  Figure  4.2c  shows  the  signal  of  another  droplet  of 
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Pi  (MIN) 


Pi  (MAX) 


fa! 


f"2  •  ^  --  *  ■  •  '■  ■  *  ■*  v  ^  >  *  • 


FIGURE  4.2.  Visibility  and  Intensity  of  Two  Size  Droplets 

(a)  and  (c)  correspond  to  droplets  traveling  thru  the 
middle  of  an  undisturbed  probe  volume;  (b)  could  be  from 
a  small  droplet  traveling  thru  the  middle  of  a  disturbed 
probe  volume  or  a  large  droplet  traveling  thru  the  edge 
of  the  Gaussian  probe  volume. 


diameter  d2  but  with  correct  visibility  V2  and  corresponding  pedestal 
Ip  •  This  last  one  will  also  be  accepted. 

Since  the  intensity  limits  are  functions  of  droplet  size,  a  vari¬ 
able  filter  like  the  one  shown  on  Figure  4.3  is  necessary.  It  contains 
the  intensity  limits  Ip(min),  Ip(raax)  corresponding  to  any  visibility 
V.  Figures  4.4a  and  4.4b  show  the  raw  data  and  software  filtered  data 
in  an  intensity/visibility  plot.  Figure  4.4a  shows  that  the  intensity 
is  bounded  by  Ip(max)  (with  a  few  exceptions  due  to  light  re-radiated  by 
upstream  droplets).  Points  with  different  intensity/visibility  combina¬ 
tions  can  be  either  due  to  droplets  moving  through  the  tail  of  the 
Gaussian  Intensity  profile  of  the  probe  volume  or  to  droplets  with  an 
erroneous  visibility.  By  limiting  the  intensity  of  the  pedestal  for 
which  a  signal  is  considered  valid,  two  things  are  accomplished:  (1) 
the  limits  of  the  probe  volume  are  carefully  established;  and,  (2)  most 
importantly,  most  of  the  signals  with  erroneous  visibility  are  rejected. 


4 .3  Results  Obtained  with  the  V/I  Technique 

Experimental  results  are  presented  to  illustrate  the  accuracy  and 
resolution  of  the  visibility/intensity  technique,  and  how  it  compares 
with  visibility  only. 


Effect  of  Beam  Blockage  on  Size  Distribution 

The  effect  of  a  very  dense  spray  blocking  the  laser  beams  before 
the  probe  volume  was  studied.  Figure  4.5a  shows  the  histogram  produced 
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whtm  mouod  1  spurned  droplets  travel  through  th,e  middle  of  an  undisturbed 
probe  volume.  Figure  4.5b  shows  a  histogram  obtained  with  the  visibil¬ 
ity  technique  of  the  same  monodispersnd  droplets  but  now  the  laser  beams 
Have  traveled  through  a  real  spfay  before  crossing*—  We  made  sure  'that 
the  spray  was  not  going  into  the  probe  volume.  The  fringe  pattern  pro¬ 
duced  by  crossing  two  laser  beams  is  well  behaved  and  theoretically  pre¬ 
dicted  in  the  undisturbed  case.  However,  when  the  beams  travel  through 
a  spray,  this  fringe  pattern  can  be  considerably  altered  due  to  random 
beam  blockage.  The  result  can  be  a  distorted  Doppler  signal  with  a  vis¬ 
ibility,  both  lower  and  higher  than  that  predicted  theoretically.  Con¬ 
siderable  error  in  the  size  distribution  will  then  be  obsorved.  The 
size  of  the  error  depends  on  the  visibility  itsslf.  It  is  quits  appar¬ 
ent  from  the  histograms  that  what  should  be  a  monodispersed  distribution 
Is  measured  as  a  very  broad  one.  The  next  histogram  shown  on  Figure 
4.5c  shows  the  results  obtained  with  V/I  when  the  beams  travel  through 
the  spray  before  crossing  (same  conditions  as  Figure  4,5b)»  Notice  that 
a  very  narrow  distribution  Is  thus  obtained. 

V/I  Measurements  with  the  Berglund-Liu  Droplet  Generator 

The  calibration  and  verification  source  used  in  these  experiments 
was  the  Burglund  Liu  monodispersc  droplet  generator.  This  generator  can 
work  in  two  modes:  (a)  it  can  generate  a  string  of  droplets  of  known 
size  and  equally  spaced  at  the  measurement  point;  (b)  using  a  cap  and 
dispersion  air,  these  droplets  can  be  dispersed  into  a  spray  of  monodis- 
peracd  droplets. 
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1c  has  been  observed  in  previous  studies^  that  under  some  disper¬ 
sion  sir  conditions,  some  of  the  droplets  wili-coLl i do  and  form  other, 
droplets  with  double  or  even  triple  the  volume.  Therefore,  if  the  pri¬ 
mary  droplet  diameter  is  d0,  then  droplets  with  diameter  2 ^d0  and  3^ 
dQ  can  be  produced.  It  must  be  pointed  out  that  in  these  studies,  we 
did  not  obtain  photographic  measurements  to  verify  the  presence  of  the 
secondary  droplets.  The  size  of  the  primary  droplets  was  accurately 
predicted  from  the  flow  rate  and  frequency  of  vibration  of  the  pinhole. 

One  limitation  imposed  by  the  Berglund  Liu  is  that  it  cannot  pro¬ 
duce  droplets  of  a  given  size  at  arbitrary  spacing.  Since  the  measure¬ 
ment  instrument  requires  the  presonce  of  only  one  droplet  in  the  probe 
volume,  that  limits  the  size  droplet  that  can  be  moasured. 

The  procedure  used  in  these  experiments  was  to  produce  a  string 
of  large  monodispersed  droplets  to  calibrate  the  instrument.  Than 
smaller  droplets  were  produced  by  increasing  the  frequency  of  vibration 
of  the  orifice,  and  with  the  dispersion  air  a  spray  of  those  droplets 
was  formed.  Typical  number  density  at  the  probe  volume  was  500  drop- 
lets/cc,  and  the  diameter  of  the  spray  at  the  plana  of  measurements  was 
about  4  mm.  ,  . . . 

A  081  system  with  V/I  capability  was  used  to  obtain  the  data" 
shown  below. 

Optical  conf igurationj 

•  Collection  angle  of  30* 

•  Transmitting  lens  ■  495  mm 

•  Waist  diameter  ■  300  pm 

•  Fringe  spacing  *  13.6  pm 

•  Collection  F if  ■  5 
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The  high  voltago  to  the  photomultiplier  was  established  with  a 
droplet  string  of  diameter  100  pm. 

The  uLxq  range  of  interest  was  9  to  92  pm. 

A  spray  of  droplets  with  primary  diameter  of  53  pm  was  then  pro¬ 
duced  using  the  Bcrglund  Liu  with  a  flow  rata  of  0.21  cc/min  and  fre¬ 
quency  of  45. 8K  Hz  and  the  dispersion  air.  For  the  conditions  selected 
above,  this  size  will  produce  a  visibility  of  58%. 

Notice  that  for  primary  droplets  of  53  pm,  the  droplets  will  have 
n  diameter  of  2^  x  53  -  66  pm,  and  the  triplets’  diameter  will  be 
31/3  x  53  ■  76  pm. 

Results  are  shown  for  both  visibility/intensity  and  visibility 
only  to  provide  some  comparison. 

Figure  4.6  shows  the  measurements  of  the  monodispersod  spray  of 

53  pm. 

Figure  4.7a  shows  the  data  of  the  spray  formed  of  primary 

droplets  and  doublets. 

Figure  4,7b  shows  results  similar  to  Figure  4,7a  but  obtained 
with  visibility  only.  Notice  how  much  broader  the  distributions  are. 

Figure  4.7c  shows  a  spray  containing  primary  droplets,  doublets 
and  triplets.  . . .  .  .  . .  ............ 

Figu  re  4,7d  shows  result'srsimll^ire'b'  but  obtained7  • 

with  visibility  only.  Notice  that  the  distribution  is  broader  end  the 
resolution  is  not  as  fine. 

Both  the  accuracy  and  resolution  of  these  measurements  are  very 
good.  The  theoretically  predicted  sizes  are  53  pm  ,  66  pm,  and  76  pm. 

The  corresponding  ineasurad  diameters  are  (52  to  54),  (63  to  65),  and 
73  pm. 
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DIAMETER  (Mm) 


VELOCITY  (m/s) 


FIGURE  4.6.  VISIBILITY/INTENSITY  MEASUREMENTS  OF  MONODISPERSED  SPRAY 
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FIGURE  4.7.  SIZE  DISTRIBUTIONS  OF  BIMODAL  AND  TRIMODAL  SPRAYS 
WITH  VISIBILITY  AND  VISIBILITY/INTENSITY 
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4 .4  The  V/I  Probe  Volume 

An  important  aspect  of  any  optical  technique  is  the  recognition 
and  determination  of  the  probe  volume.  The  implication  is  that  differ¬ 
ent  size  particles  scatter  different  amounts  of  light  (typically  propor¬ 
tional  to  the  square  of  the  diameter).  Also,  the  intensity  of  the  probe 
volume  is  non-uniform  (typically  Gaussian).  Therefore,  for  a  given 
threshold  level  the  volume  of  detectability  is  a  function  of  particle 
size.  This  will  introduce  a  natural  bias  in  the  histogram  of  particle 
size.  To  correct  this  bias  the  probe  volume  needs  to  be  computed  as  a 
function  of  particle  size  arid  optical  configuration  and  then  used  to 
correct  the  biased  histogram. 


A 


CF 

D 


P 

Dlaser 


area  of  measuring  region 
waist  radius 
collection  factor 
diameter  of  pinhole 
diameter  of  laser 
diameter  of  droplet 

maximum  droplet  diameter  (dQ/6CF)  =  .845)  for  a  given  size 
range 

focal  length  of  transmitting  lens 

focal  length  of  1st  collecting  lens 

focal  length  of  2nd  collecting  lens 

ground  to  peak  threshold  (front  panel  in  volts) 

dc  saturation  (0.65  volts) 

fringes  required  by  processor  (front  panel) 

pedestal  intensity 

visibility  of  size  d 

visibility  of  size  dQ  (VQ  ■  .151) 

coordinate  normal  to  fringes 


y,z 

6 


coordinates 
fringe  spacing 
Y  :  crossover  angle 

9  :  angle  of  collection 

4 .5  Simplified  Probe  Volume  Model 

The  cross-sectional  area  is  given  by: 


sin9  f^ 


where  y  must  be  calculated  based  on  the  Gaussian  profile  and  d  (droplet 
diameter) . 

Let  us  call  Ig  the  saturation  level  Imposed  by  the  electronics. 
Then  for  the  largest  detectable  droplet: 

Ipeak  <  Ig. 

The  intensity  profile  of  a  probe  volume  formed  by  crossing  two  Gaussian 

beams  is  given  by: 


2  2 


I  -  21  exp  {  (-—)  [x2+y2  +  } 


„  4*  x  sinf^-) 

r  .  x  zYi  . 

[cos  h  [ - » — LJ  +  cos  - \ -  J 

b 

o 
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For  a  pinhole  limited  signal,  we  can  assume  z  to  be  small;  therefore. 


I  -  2Iq  exp  [  (-^r)  (x2+u2)  ]  •  [1  +  cos  ] 

b 

o 


The  pedestal  scattered  by  the  droplet  of  diameter  d  Is  given  by: 


P  "  2IoKod2  exp  t  ^  ^*2+y2^  1 

b 

o 


The  intensity  scattered  by  the  same  droplet  is  given  by: 


21,  V2  exp  t  (^T)  t*2^2)  )  t  1  +  SSSTrai  •  v  ) 


(2) 


The  peak  intensity  occurs  at  x*0,  y*0;  and  for  the  largest  droplet, it  is 
given  by: 


Ipeak  -  21  K  d  2(l  +  V  )  <  I 
r  O  O  O  v  O'  8 


Solving  for  1Q: 
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(3) 


I 

o 


I 

s 


■2K  d  2  (l  +  V  ) 
o  o  v  o' 


The  expression  for  the  peak-to-peak.  i-ntens i  ty  I^”:;can  be  obtained '  by'  .  . 
setting  the  cosine  term  to  +1  and  -1  in  Equation  2,  and  then 
subtracting: 


Substituting  Equation  3  in  the  above  equation,  we  get: 


(4) 


and 
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In  the  old  probe  volume  model,  it  was  .sufficient  that  I  >  I  .  . 

ac  min 

The  new  intensity/visibility  software  filter,  however,  limits  the  inten¬ 
sity  further  by  limiting  the  pedestal  of  acceptability. 

-  '■  — T'.'u  .  ..  .  --  -----  •’  ••  '.-...•a--*’-***  wprn  u&mrj-Mi 

To  criterion  that  I  >  Im^n  is  absolutely  necessary  but  not  suf¬ 
ficient,  except  for  the  first  A  bins  which  have  no  minimum  pedestal 
intensity  limitation. 

The  intensity/visibility  software  filter,  which  is  contained  in 

2 

Table  1,  was  constructed  such  that  is  proportional  to  d  . 

’  max  r  v  max 

For  instance  for  bin  35,  d_av  ■  .595  and  P,_a„  ■  271.  For  bin  8, 

’  max  max 


max 


.2185. 


P 

maXg 


r  .2185>2 
' .595  > 


x  271 


36.5 


Given  that  the  visibility  falls  in  any  one  bin,  the  actual  pedestal 
scattered  by  the  droplet  has  to  fall  between  and  Praax>  where  PMX 

corresponds  to  the  droplet  moving  through  the  center  x*y»z*0of 
the  probe  volume.  Given  this  visibility  V  there  is  a  bin  //  associated 
which  corresponds  to  a  diameter  d  (as  given  by  the  visibility  Table  2). 
The  intensity/visibility  software  filter  (included  in  Table  1)  imposes 
the  condition: 
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TABLE 

1. VISIBILITY 

INTENSITY  TABLE 

. - 

^  ...  .  . 

V". 

IMIN 

IMAX 

XMIN 

*MAX 

BIN  NO. 

(mV) 

(mV) 

BIN  NO. 

(mV) 

(mV) 

1 

0 

16 

40 

225 

348 

2 

0 

17 

41 

235 

362 

3 

0 

20 

42 

247 

378 

4 

0 

23 

43 

258 

394 

5 

4 

25 

44 

269 

410 

6 

6 

27 

45 

280 

427 

7 

9 

33 

46 

293 

443 

8 

11 

37 

47 

305 

460 

9 

14 

42 

48 

317 

479 

10 

18 

46 

49 

329 

495 

11 

21 

51 

50 

343 

515 

12 

25 

57 

51 

355 

533 

13 

28 

63 

52 

368 

550 

14 

33 

69 

53 

382 

571 

15 

37 

75 

16 

42 

82 

17 

47 

90 

18 

52 

97 

19 

57 

105 

20 

63 

114 

21 

69 

122 

22 

75 

131 

23 

81 

139 

24 

87 

149 

25 

95 

160 

26 

101 

170 

27 

108 

180 

28 

116 

191 

29 

124 

203 

30 

132 

214 

31 

140 

225  . 

32' 

—  149  •.  T*:--..- 

237 '  — 

33 

158 

251 

34 

167 

264 

35 

175 

277 

36 

184 

290 

37 

194 

304 

38 

205 

318 

39 

215 

332 
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TABLE  2.  THE  VISIBILITY  FUNCTION 


BIN 

vis 

VIS 

REL 

DIA/ 

BIN 

VIS 

VIS 

REL 

DIA/ 

NO. 

LOW 

HI 

BW 

DELTA*CF 

NO. 

LOW 

HI 

BW 

DELTA*CF 

1 

.986 

.989 

4 

.0813 

40 

.385 

.404 

20 

.6546 

2 

.981 

.985 

5 

.096 

41 

.365 

.384 

20 

.6692 

3 

.975 

.98 

6 

.1107 

42 

.345 

.364 

20 

.684 

4 

.968 

.974 

7 

.1254 

43 

.326 

.344 

19 

.6987 

5 

.961 

.967 

7 

.1401 

44 

.306 

.325 

20 

.7134 

6 

.953 

.96 

8 

.1548 

45 

.287 

.305 

19 

.728 

7 

.944 

.952 

9 

.1695 

46 

.268 

.286 

19 

.7428 

8 

.935 

.943 

9 

.1842 

47 

.249 

.267 

19 

.7575 

9 

.924 

.934 

11 

.1989 

48 

.231 

.248 

18 

.7722 

10 

.914 

.923 

10 

.2136 

49 

.213 

.23 

18 

.7868 

11 

.902 

.913 

12 

.2283 

50 

.195 

.212 

18 

.8016 

12 

.89 

.901 

12 

.243 

51 

.177 

.194 

18 

.8163 

13 

.877 

.889 

13 

.2577 

52 

.16 

.176 

17 

.831 

14 

.864 

.876 

13 

.2724 

53 

.143 

.159 

17 

.8456 

15 

.85 

.863 

14 

.2871 

16 

.835 

.849 

15 

.3018 

17 

.82 

.834 

15 

.3165 

18 

.804 

.819 

16 

.3312 

19 

.788 

.803 

16 

.3459 

20 

.772 

.787 

16 

.3606 

21 

.755 

.771 

17 

.3753 

22 

.737 

.754 

18 

.39 

23 

.72 

.736 

17 

.4047 

24. 

.701 

.719 

19 

.4194 

25 

.683 

.7 

18 

.4341 

26 

.664 

.682 

19 

.4488 

27 

.645 

.663 

19 

.4635 

28 

.626 

.644 

19 

.4782 

29 

-.606 

.625 

20 

.4929  . 

30 

.587 

.605 

19 

.5076 

31 

.567 

.586 

20 

.5223 

32 

.547 

.566 

20 

.537 

33 

.527 

.546 

20 

.5516 

34 

.506 

.526 

21 

.5664 

35 

.486 

.505 

20 

.5811 

36 

.466 

.485 

20 

.5958 

37 

.446 

.465 

20 

.6104 

38 

.425 

.445 

21 

.6252 

39 

.405 

.424 

20 

.6399 
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where  P  is  measured  by  a  peak  detector  of  the,  Gaussian  envelope  and, 

therefore,  it  corresponds  to  x*0.  I  and  P„v  correspond  to  each 

s  I  ma  x 

other. 


Given  a  bin  //,  there  is 
a  correspondence  between 
diameter  and  pedestal. 


drain 

d 

max 

P  ,  -  0.7  Pfd  .  )/ 

min  v  rain'  max 

P 

max 

where: 


min'  max 


drain  81  ® 

— - - e  , 

do  ^1+Vo) 


p 

rain 


0.7  d  .  I 

min  s , 


V  O+vJ 


-  ■  ■■  d  2(l+V  )p  ,  d  2(  1H  l.P  .. 

Therefore,  I  -  "  —  °'-  fn  =  ° 

S1  0.7  d  ,  2  d  2 


f&i-  >*»*'** 


min 


max 


Substituting  for  I  above: 

®  1 


d2  P 

„  max 

P  2  exp  v  2  >  "  ‘min 

a  d 

max  o 


> 


Solving  for  y 
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2 


b  p  , 

—  ln  l VJ J 

Kd 


where  K  = 


max 


max 


at  any  bln  # 


(6) 


For  instance  for  bin  #35,  K 


.271 


.595 


.765 


Note  that  the  value  of  y  will  only  remain  constant  if  the  characteristic 


size  is  d  .  ,  in  which  case 
min 


However,  I  have  chosen  dcharact  to  be  the  value  given  in  the  visibility 
table,  which  is  larger  than  d[n^n;  therefore  y  >  0.42.  For  the  high  vis¬ 
ibility  bins,  dcharact  »  ^min’  and>  therefore,  we  should  expect 
y  »  .42  as  shown  in  the  example.  The  expression  for  y  based  on  the 
threshold  and  saturation  intensity  is  given  by: 
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For  the  conditions: 


ITh  =  .0005V 

Is  -  .63V 

<*o  =  .845 

VQ  =  .151 


Also  from  Table  1  at  bln  #  35,  d'  =  .595  and  Pmax  =  .271V  at  y  =  0. 


.27  if  .845) 2  fl.151)  _  >63v 
.5952 


K 


.271 


.595 


.765 


Therefore,  Equation  6  yields: 
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where  Pm^n  is  given  in  volts 


Equation  7  yields: 


y  <  b 

o 


6.3  x  10 

'2 

d  V 


-4 


1 


For  example: 


-• 


d' 

V 

P  Zl6-) 

min  b 

o 

y(7)/bo 

53 

.8456 

.15 

.373  .436 

1.25 

‘  • 

35 

.5811 

.49 

.171  .44 

1.34 

13 

.2577 

.88 

.0278  .55 

1.12 

5 

.1405 

.965 

.0042  .8 

.84 

4 

.1254 

.97 

0  0 

.77 

Notice  that 

the  relationship 

between  y  and  bin  //  would 

be  a  constant  if 

*• 

for  Pmi.  n  I 
min 

used  the 

corresponding  dm^n  instead  of  the 

d  given  on  the  V 

Table. 


_• 
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5.0 


THE  IMAX  TECHNIQUE 


The  method  discussed  here  bases  the  size  measurement  on  the  abso¬ 
lute  intensity  scattered  by  the  droplet  crossing  the  probe  volume,  and 
the  velocity  measurement  on  the  classical  Doppler  signal.  It  is 
referred  to  as  the  IMAX  method.  In  situ  single  particle  counters  are 
limited  because  of  the  nonuniforra  profile  (typically  Gaussian)  of  laser 
beams.  Under  this  condition  a  particle  crossing  the  middle  of  the  beam 
will  scatter  more  light  than  a  similar  particle  crossing  through  the 
edge.  Therefore,  the  relationship  between  size  and  scattered  light  is 
not  unique. 

To  circumvent  this  problem,  two  beams  of  unequal  size  are  crossed 
such  that  the  small  beam  identifies  the  middle  of  the  large  beam  and 
therefore,  removing  the  Gaussian  ambiguity.  Laser  beams  are  chosen 
because  of  their  spatial  and  temporal  coherence  and  because  of  their 
size.  These  beams  will  interfere  where  they  cross  and  a  fringe  pattern 
will  be  formed  in  the  middle  of  the  large  beam.  Signals  exhibiting  an 
ac  modulation  will  have  crossed  the  fringe  pattern  and  therefore,  the 
middle  of  the  large  beam.  Both  size  and  velocity  of  individual  spher¬ 
ical  particles  can  be  extracted  from  this  signal. 

At  least  two  approaches  can  be  used  to  implement  the  above  con¬ 
cept.  The  first  consists  of  crossing  two  laser  beams  of  different  diam¬ 
eters  but  with  the  same  wavelength.  In  the  second  two  small  beams  of 
one  wavelength  cross  in  the  middle  of  a  larger  beam  of  another  wave¬ 
length.  These  approaches  will  now  be  explained  starting  with  the 
second . 
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Figure  5.1  illustrates  the  probe  volume  of  this  method.  A  spher¬ 
ical  particle  crossing  through  the  fringes  will  also  cross  through  a 
region  of  almost  uniform  intensity  of  the  large  beam.  The  ratio  of  the 
two  beam  diameters  will  establish  the  uniformity  of  the  intensity  inci¬ 
dent  on  the  droplet. 

If  we  refer  to  the  small  beam  as  1  and  the  large  beam  as  2,  the 
intensity  profiles  in  the  probe  volume  can  be  spectrally  separated,  and 
given  by: 


I,  =  21  exp  {(-5~)[x2+y2+z2Y2/4] 


(cosh  (^) 


+  cos 


4rrxsin(^-) 


].  (8) 


and 


l2  =  lo  exp  t  “  *T~^x2+y2^  '  (9) 

2  b 

°2 

Where  IQ  is  the  center  intensity,  y  is  the  intersection  angle,  bQ  the 
waist  radius,  A  the  laser  wavelength,  and  x,y,z  the  coordinates.  The  z 
dependence  of  the  large  beam  is  neglibible.  If  we  also  assume  that 

jV 

2~  »  0(which  is  an  excellent  assumption  since  a  pinhole  in  the  receiver 
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will  limit  the  value  of  z),  the  intensity  scattered  by  a  spherical  par¬ 
ticle  is  given  by: 


I 

S1 


21  K 
o .  o 


d2  exp  +  COS  2^ 

b  x 

°1 


V]  , 


(10) 


and 


I 

S2 


I  K  d2 
°2  °2 


exp 


(ID 


where 


Ko 


?/. 


'[e  (6.n)  0(6)]  +  [E2(0,n)  D(0)] 

refraction  reflection 


(asinG) 


dA 


sin  9  diff . 


(12) 


Kq  is  a  scattering  coefficient,  D  is  the  divergence,  n  the  index  of 
refraction,  e  the  fraction  of  energy  for  every  ray  reflected  or 
refracted,  d  is  the  diameter  of  a  spherical  particle,  V  is  its 
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visibility,  r  the  distance  from  the  scattering  center  to  the  collecting 


lens,  and  9  the  collection  angle  measured  from  the  forward  direction. 

It  can  be  shown^^  that  the  contributions  of  these  three  terms  is  a 
function  of  angle.  To  keep  the  following  analysis  simple  we  will  choose 
a  collection  angle  9  of  30°.  Assuming  that  the  solid  angle  of  the 
collecting  lens  is  small  with  respect  to  0  we  can  then  compare  the 
intensities  scattered  at  the  discrete  angle  0. 

The  refraction  or  reflection  terms  are  given  by: 

2  2„ 

iR  =  a  e  D  , 

and  the  diffraction  by: 


~~  (ctsin0) 


sin  0 


where  a  -  y—  is  the  size  parameter,  and  Jj  the  Bessel's  function  of  the 
first  kind.  Reference  9  shows  that  for  S  polarization: 


i(refraction)  =  1.0375 
l(ref lection)  =  0.0785 


and  it  can  be  calculated  that 
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ipjCa)  <■  .08^  a" 


for  a  >49  (d  >  7.6  pm).  At  this  angle  the  combined 


reflected  ana  diffracted  light  represent  16£  of  the  refracted  light. 

Which  implies  that  in  general  the  three  components  must  be  taken  into 

account.  It  is  assumed  that  the  collecting  lens  is  in  the  plane  of 

symmetry  of  the  crossing  beams.  The  above  equations  show  that  the 

scattered  light  intensity  is  proportional  to  the  square  of  the  diameter 

of  the  sphere.  This  approximation  is  only  valid  for  diameters  larger 

than  about  10  wavelengths.  Note  that  the  angle  8  (deviation  of  incident 

pencil  of  light)  to  an  arbitrary  position  of  the  lens  is  different  for 

both  beams  (except  on  the  line  of  symmetry).  However,  the  integrated 

value  is  the  same  as  long  as  the  collecting  lens  is  in  the  plane  of 

svmmetry.  In  this  technique  the  ac  modulation  of  1  is  used  to  estab- 

S1 

lish  detectability  of  the  signal  and  to  measure  the  velocity.  For  a 
given  threshold  level,  the  larger  the  dynamic  size  range,  the  farther 
from  the  center  of  the  beam  the  largest  particle  can  cross  and  still 
exhibit  ac  modulation  and,  therefore,  be  detectable.  This  will  intro¬ 
duce  an  error  in  the  measurement  as  will  be  discussed  in  the  next 
section. 

The  peak  to  peak  ac  signal  is  then  given  by: 


I 

ac 


1 


41  K  Vd2exp  [(-  ~-)(xV)] 
°L  b* 


(13) 
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and  the  pedestal  is  given  by: 


P,  =  21  K 
1  °1  °1 


d2  exp  [(  -  -y-)(x2+  y2)] 
b 

°1 


Dynamic  size  range  and  error  analysis 

There  are  two  known  sources  of  error  in  the  IMAX  method.  First 
is  the  error  introduced  by  beam  blockage  produced  by  other  particles  in 
the  trajectory  of  the  beams  and  scattered  light.  This  one  is  difficult 
to  quantify  and  will  be  discussed  with  the  results.  Second,  since  the 
small  beam  has  a  finite  diameter  the  particles  can  travel  a  small  dis¬ 
tance  away  from  the  center  of  the  big  beam  and  still  cross  sufficient 
fringes  to  be  detectable.  The  farther  from  the  center  it  can  be  detect¬ 
able,  the  larger  the  error  it  can  produce.  It  will  be  shown  here  that 
this  error  increases  with  dynamic  size  range,  and  decreases  with  ratio 
of  the  large  to  small  beam. 

The  dynamic  range  is  limited  by  the  electronic  and  optical  noise. 
In  general,  we  can  establish  that  there  is  a  threshold  to  the  minimum 
processible  ac  signal,  and  a  saturation  level  associated  with  the 
largest  signal  processible  by  the  electronics.  All  the  processible  sig¬ 
nals  must  fall  between  these  two  limits.  Equation  13  gives  the  peak  to 
peak  ac  signal.  Notice  that  both  the  diameter  (d)  and  the  visibility 
(V)  influence  the  above  expression. 

To  simplify  this  analysis  we  will  neglect  the  contribution  of 
diffraction,  that  is  the  collection  angle  will  be  off-axis.  Single 
particle  counters  usually  collect  off-axis  light  (6  >  30°)  for  the  sake 
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of  producing  a  small  probe  volume.  Therefore,  the  above  assumption  is 
very  reasonable  for  spheres  larger  than  7  pm.  The  dynamic  range  is  act¬ 
ually  larger  when  diffraction  is  included. 

The  minimum  processlble  ac  signal  is  produced  by  the  smallest 
particle  (with  diameter  d^^  crossing  the  probe  volume  at  y  *  o  and 
crossing  the  minimum  number  of  fringes  required  by  the  electronics 
(x  =  bDi). 

Therefore,  lac  =  4l0^K0^dmin  exP  ("2)  where  it  has  been 
assumed  that  the  visibility  of  the  smallest  particle  is  1. 

The  maximum  ac  signal  is  produced  by  a  droplet  of  diameter  dQ 
(not  necessarily  d„  1  and  visibility  V  . 

tuaX  U 


‘ac  =  4Io,Vo2vo  6XP  ■J“)(x2+y2^]  * 

max  11  b 

°i 


(15) 


The  farthest  from  the  center  this  particle  can  be  detected  is  when: 


^ac,„,„  "  ‘ac  .  „ 
max  mm 


and  it  crosses  sufficient  number  of  fringe  (x  =  bo^).  Therefore, 


41  K  d  2  exp  [( - )(b2  +y2)]  =41  K  d2  exp  (-2)  . 

Oj  ol  o  K  LV  b2  Oj  J  Ji  Oj  Oj  min 

°1 


-40- 


84-2286-06/36 


Solving  for  the  size  range,  we  obtain: 


d 

= 

'  d  .  1 
min 


_  exp(-2) 
V 


and 


2(b^  +y2) 

exp  [ - - - ]  = 


'1 


1  (hi) 

—  exp 


,  2  2  „  rr  o  ^2 

2y  =  b  tnlCj—) 

I  min 


- 7~]  -2  b2 

exp(-2)J  o | 


which  results  in 

2y2  -  bo  *«[v0(jr-)2 

l  min 


(16) 


The  error  in  the  size  will  be  due  to  an  error  in  the  intensity  scattered 
by  the  large  beam  when  the  particle  crosses  at  y  >  o  (notice  that  the 
particle  must  cross  through  x  =  0).  If  we  state  that  the  large  beam  is 
m  times  larger  than  the  small  one, 


To  simplify  the  error  analysis,  let  us  assume  that  the  optical  param¬ 
eters  are  chosen  such  that  for  any  size  range  the  visibility  of  the 
largest  droplet  is  larger  than  zero. 

That  is  V(d  ,)  >  0  and  d'  <1  where  d1  =1  produces  the 
max  max  max 

first  zero  in  the  visibility.  Then  the  maximum  ac  signal  is  produced  by 
a  particle  of  diameter  dQ  -  0.626  dmax  and  its  visibility  is  43%. 
Substituting  these  values  into  Equation  16,  we  obtain: 


— ^ —  =  Jin  [.1685  ( 
b 

°1 


d  o 
max\  2 

d  .  ‘ 
min 


The  maximum  errors  resulting  from  the  finite  size  of  the  beams  can  now 
be  presented  in  tabular  form. 


m 

^raax^rain 

2 i_ 
b2 

1 

S2 

I  K2 

max  error 

°i 

°2  °2 

7 

10 

2.82 

.94 

32 

7 

20 

4.21 

.92 

4% 

7 

30 

5.02 

.90 

5% 

5 

10 

2.82 

.89 

5% 

5 

20 

4.21 

.85 

87. 

5 

30 

5.02 

.82 

10% 

5.2  Single  Color  System 

Figure  5.2  shows  a  schematic  representation  of  the  probe  volume 
formed  by  crossing  two  beams  of  the  same  wavelength  but  different 
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diameters.  There  are  several  fundamental  differences  between  this  and 
the  two  color  system:  1)  Fringes  are  formed  everywhere  the  large  and 
small  beams  mix;  2)  The  pedestal  of  the  large  beams  can  not  be  spec¬ 
trally  separated  but  electronically  filtered.  The  intensity  distibution 
of  either  beam  can  be  expressed  as: 


1  =  I  exp  [-  ~~  (x2+y2)] 

1  b 


(18) 


and 


l2  =  ro  exP  [  — |—  C^2+y2)]  •  (19) 

2  b 

°2 

The  intensity  scattered  by  the  spherical  particle  can  be  expressed  by: 


1=1  +1  +2/1  I  cosB’V 

S  S1  S2  V  S[  ®2 


where 


I  =  IK  d2 
s  o 
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FIGURE  5,2.  PROBE  VOLUME  OF  SINGLE  COLOR  IMAX  TECHNIQUE. 


8  is  the  phase  angle,  and  again  it  is  assumed  that  the  collecting  lens 
is  on  the  plane  of  symmetry  of  the  two  crossing  beams.  The  pedestal  is 
given  by 


P  =  I  +  I 
S1  s2 


and  the  peak  to  peak  ac  intensity: 


I 

ac 


=  2  fi  l  V 

V  S1  S2 


(cos0  =  1)  -  2 Jl  I  V 


(cos0  =  -l) 


=  4  fi  r 

V  S1  £ 


That  is. 


-f“  (x2+y2)]  +  l0  K0d2  exp  [  -  ~  (x2+y2)] 
b  °2  0  b 

°1  °2 


and 


I 

ac 


4K  R  I  d2  V  exp  f 
of  o{  o2 


/  2  2  ■> 

_  l.x.  JX  J, 


(21) 


(20) 
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Following  che  same  approach  as  before,  and  again  neglecting  the  contri¬ 
bution  of  diffraction  for  the  off-axis  collection  angle.  If  we  let 


m  *  /b.  ,  we  obtain: 

°2  °1 


I  =  4K  fl  T~  d  exp  [  -1  -  -4] 

ar  r\  /  mi  n  L  /  J  ' 


ac  .  o  /  o„  o„  min 

min  J  22 


I  =  4K  /i I  d  V  exp  f 

ac  o  /  o,  o_  o  o  L 

max  V  1  2 


fbol+y2) 


(b^V) 


The  variation  in  y  will  result  in  an  error  of  the  measured  pedestal.  If 
we  can  assume  that  m  is  large,  then  the  pedestal  of  the  large  beam  can 
be  recovered  by  electronic  filtering  and  it  is  given  by  the  second  term 


of  Equation  20.  This  variation  in  y  is  obtained  equating  22  and  23,  and 
we  obtain 


-  -  *•  • 


2k2  „ 

m  b  2 

~  o.  d  V 

2  1  .  r  o  oi 

y  =  —p - In  L — 5 — j 

m  +1  d  7 

min 


The  error  in  the  size  will  result  from  droplets  crossing  the  pedestal  of 
the  large  beam  at  x  =  o  and  y  given  by  24.  That  is: 


I  K  d 


[-rH 

mb 

°1 
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Making  similar  assumptions  as  before: 


d  =  .626  d  and  V  =  0.43 
o  max  o 


The  maximum  errors  resulting  from  the  finite  size  of  the  beams  can  be  pre¬ 
sented  in  tabular  form: 


ra 

dnax^dmin 

_2)L_ 
2.2 
m  b 

°1 

1 

S2 

I  K  d2 
°2  °2 

max  error 

7 

10 

.113 

.89 

5% 

7 

20 

.168 

.84 

8% 

7 

30 

.2 

.82 

10% 

5 

10 

.22 

.80 

10% 

5 

20 

.32 

.72 

15% 

5 

30 

.39 

.68 

18% 

It  should 

be  pointed  out 

that  the 

technique  discussed  here  bases  the 

detectability  of  the  sig 

nal  on  the 

modulated  or 

ac  part  of  the  scattered 

light.  It  is,  therefore,  important  that  no  zeros  are  present  in  the 
visibility  function. 
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SPC's  based  on  absolute  scattered  light  require  calibration  at 
one  datum  point.  When  measuring  dense  sprays  beam  attenuation  will  mat'' 
this  calibration  difficult.  An  algorithm  is  discused  here  that  will 
eliminate  this  problem. 

The  algorithm  is  based  on  measuring  the  visibility  of  any  size 

/ 1 6  \ 

class  droplet.  Combining  the  visibility  and  intensity'  '  of  the  signal 
a  very  accurate  measurement  of  such  size  can  be  made.  Knowing  the  size 
of  the  chosen  corresponding  visibility  one  can  then  establish  the  ampli¬ 
tude  of  the  pedestal  of  this  signal  by  adjusting  the  gain  to  the  photo 
detector.  The  logic  is  as  follows:  for  a  given  visibility  (chosen  in 
the  most  accurate  region  of  the  visibility  function)  there  is  an  asso¬ 
ciated  size  and,  therefore,  an  intensity  scattered  from  the  large  beam. 
Since  the  intensity  of  this  large  beam  is  almost  uniform  the  scattered 
light  associated  with  the  chosen  visibility  will  be  almost  constant  in 
the  absence  of  errors  in  the  visibility.  There  are,  however,  errors 

associated  with  measuring  visibility  in  a  dense  spray  and  these  points 
(16) 

are  rejected  by  V/Iv  ’ .  The  process  is  an  iterative  one  where  V/I  is 
implemented  by  choosing  several  narrow  intensity  bands  to  establish 
points  out  of  control.  The  gain  to  the  photo-detector  is  adjusted  until 
the  majority  of  the  points  fall  in  the  pre-established  intensity  band 
corresponding  to  the  measured  size. 

5. A  The  IMAX  Probe  Volume 

The  dynamic  range  exercise  indicated  that  droplets  that  scatter 
light  with  large  modulation  are  detectable  over  a  larger  region  than 


-48- 


84-2286-06/36 


those  that  scatter  light  with  little  modulation.  As  a  result,  the  probe 
volume  or  region  of  detectability  is  a  function  of  droplet  size. 
Algorithms  have  been  presented^ ^  to  define  the  probe  volume  in  LDV 
type  systems.  A  very  simple  algorithm  will  be  presented  here  for  the 
two  color  system  and  off-axis  collection. 

Notice  that  only  the  small  beams  define  the  probe  volume  in  this 
case  and,  therefore,  the  probe  volume  will  be  a  standard  LDV  type 
volume.  The  approach  followed  will  be  very  similar  to  the  one  described 
in  Reference  (17). 

There  are  two  signal  levels  imposed  by  the  electronics:  the 

saturation  (I  )  and  the  threshold  (I  .  ).  The  intensity  to  the  peak 
ma  x  in  1  n 

of  the  ac  modulation  (Ipea^)  of  the  largest  signal  must  be  less  than 
I  ,  and  the  smallest  detectable  signal  must  be  larger  than  Imin 
(defined  peak  to  peak). 


I  ,  =  21  K  d2  (l+V  )  <  I  .  (27) 

peakj  0^  Oj  o  '  o'  max 

where  dQ  produces  the  largest  signal  within  the  detectable  size  range, 
and  Vq  is  its  corresponding  visibility. 

Solving  for  IQ  we  obtain: 


I 


I 

max 


2K  d 
0^  c 


( i+vo) 


(28) 
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This  expression  can  now  be  substituted  in  Equation  13  to  obtain: 


2d  V  I  _  .  „ 

,  =  — — — “  exp  [(  -  —  )  (x  +y  )] 

"1  d  fl+V  )  b 

o  K  oJ  O  , 


and  the  detectability  criterion  establishes  that  I  >  I  .  and  that 

J  acj  rain 

sufficient  fringes  are  crossed.  If  we  assume  that  the  number  of  fringes 
processed  by  the  electronics  is  equal  to  that  in  the  waist  diameter, 


x  =  b 

°1 


Solving  for  y  in  Equation  29  we  obtain: 


V  rr*in 

2  £n  21 


d  (l+V  ) 
f_0\  o  1 

V  ri>  u  J 


The  cross-sectional  area  of  sensitivity  for  an  off-axis  pinhole  limited 
signal  can  be  approximated  by: 


A(d)  - 


2v  ii 

sin9  f. 


where  Dp  is  the  pinhole  diameter  and  F ^ / f 2  determines  the  magnification 
of  the  receiving  optics.  The  effect  of  the  pinhole  and  the  coordinate 


system  are  illustrated  on  Figure  3.3. 
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As  seen  by  Equation  30  and  31  the  area  of  detectability  (in  the 
y,z  plane)  is  only  a  function  of  droplet  size  and  corresponding  visibil¬ 
ity.  This  area  can  now  be  normalized  and  provide  a  weighting  function 
to  size  distributions  otherwise  biased  to  the  large  signals. 

This  weighting  factor  is  given  by: 


w(d) 


A(d) 

A 

max 


The  number  density  of  any  size  droplet  can  also  be  calculated  by 
dividing  the  number  of  counts  of  the  given  size  by  its  corresponding 
sampling  volume  (SV). 

The  sampling  volume  is  a  cylinder  with  cross-section  A(d)  and 
length  equal  to  the  droplet  velocity  (U(d))  multiplied  by  the  effective 
sampling  time  (T). 

M 

That  is  SV(d)  =  A(d)U(d)T,  and  T  =  TS  -  £  N  tD^  -  Mr 

i-1  1 

where  TS  is  the  total  sampling  time  to  collect  M  samples,  N  is  the 
number  of  fringes,  tD^  is  the  Doppler  period  of  sample  i  and  xQ  is  the 
processor's  dead  time  after  every  sample. 

The  number  density  of  size  d  is  then  given  by: 


ND(d) 


N(d) 

A  U(d)T 
max 


(32) 
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r\ 


where  a(d)  is  the  raw  count  of  size  d.  Because  of  the  multivalued  char¬ 
acteristic  of  the  visibility  curve,  VQ  should  be  larger  than  the  first 
zero  of  the  visibility  function.  Equation  (10)  defines  the  abscissa  of 

this  function  as  d*  =  777=-  .  This  value  should  be  kept  smaller  than 

oCF 

one,  however  physical  constraints  render  this  not  always  possible.  The 

collection  factor  (CF)  can  be  easily  changed  almost  linearly  by  changing 

the  collection  F #.  In  the  experiments  reported  here  this  was  done  by 

placing  a  mask  on  the  path  of  light  scattered  by  the  small  beams.  As  an 

example,  the  spray  data  shown  on  Figures  5.9A  and  5.9B  correspond  to 
—3 

Imin  =  .6  x  10  V,  Imax  *  .7  V  (these  intensities  are  read  directly  in 

volts).  Since  very  little  data  was  measured  above  160  pm,  the  mask  was 

chosen  to  obtain  d'  =  1  at  about  160  pm.  Therefore 
'  200  200 

d  =  7-7=-  =  7777  =  1.24.  The  following  values  are  then  obtained  for  the 
o  our  iou 

nondimensional  diameter  (d ' ) ,  its  visibility  (V)  and  the  normal  probe 
volume  dimension  (y,  as  given  by  Equation  30). 


d 

20 

40 

60 

80 

100 

120 

140 

160 

180 

200 

d' 

.12 

.25 

.37 

.5 

.62 

.75 

.87 

1.0 

1.12 

1.24 

V 

.97 

.88 

.76 

.6 

.43 

.27 

.12 

.02 

.08 

.12 

y 

.35 

.9 

1.06 

1.14 

1.17 

1.15 

1.03 

.55 

1.05 

1.19 

Similar  calculations  were  performed  for  the  size  ranges  5  to  50  pm,  and 
10  to  100  pm. 
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Description  of  IMAX  Breadboard 

Figure  5.4  shows  a  schematic  diagram  of  the  two-color  breadboard 
system.  An  argon-ion  laser  provides  the  light  source.  The  beam  colors 
are  separated  by  the  dispersion  prism  and  two  of  these  colors  (the  5145A 
shown  by  the  broken  line,  and  the  4880A  shown  by  the  solid  line)  are 
used  to  define  the  probe  volume.  A  beam  expander  formed  by  lenses  LI 
and  L2  define  the  beam  ratio  parameter  m.  A  compensated  beamsplitter 
splits  the  blue  beam  into  two  beams  and  the  transmitting  lens  L3  focuses 
and  crosses  these  two  blue  beams  in  the  middle  of  the  green  beam.  The 
receiving  optics  use  a  dichroic  mirror  to  separate  the  two  colors  of  the 
scattered  light.  The  probe  volume  is  imaged  on  pinholes  in  front  of  two 
PMT's  and  interference  filters  are  used  to  perfectly  separate  the  two 
colors.  The  outputs  of  the  two  PMT’s  are  then  electronically  processed 
as  shown  on  Figure  5.5,  and  information  of  the  size  and  velocity  of 
individual  spheres  crossing  the  probe  volume  is  thus  obtained.  An  elec¬ 
tronic  processor  was  developed  for  this  purpose,  and  a  dedicated  micro¬ 
processor  was  interfaced  to  store,  display  and  analyze  the  acquired 
data.  Sample  rates  of  several  KHz  are  possible  with  this  system. 
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Results  Obtained  with  the  IMAX  Technique. 
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FIGURE  5,4.  SCHEMATIC  I  MAX  BREADBOARD  SYSTEM 


IMAX  Measurements  with  the  Berglund  Liu  Droplet  Generator 


A  vibrating  orifice  droplet  generator  was  used  to  produce  strings 
and  sprays  of  known  size  droplets.  The  procedure  used  in  these  experi¬ 
ments  was  to  produce  a  string  of  large  raonodispersed  droplets  to  cali¬ 
brate  the  instrument.  Smaller  droplets  were  then  produced  by  increasing 
the  frequency  of  vibration  of  the  orifice,  and  with  the  dispersion  air  a 
spray  of  these  droplets  was  formed.  The  spray  angle  was  about  10°  and 

O 

the  number  density  was  typically  500/cm  .  The  optical  configuration 
consisted  of  a  transmitting  lens  of  711  mm,  m  =  7,  blue  waist  diameter 
of  100  pm,  receiving  angle  of  20°,  and  receiving  lenses  of  300/495  mm. 
The  calibration  point  was  provided  by  a  string  of  110  pm  droplets  pro- 

O 

duced  with  a  flow  rate  of  0.21  cm  /min  and  a  frequency  of  5  KHz.  Figure 
5.6  shows  the  measured  size  and  velocity.  A  spray  of  primary  droplets 
of  49  pm  (produced  with  a  flow  rate  of  0.21  cm  /min  and  frequency  of 
56.9  KHz)  was  then  produced  with  the  dispersion  air.  Figure  5.7a  shows 
the  measurements  of  the  spray  of  primary  droplets.  Figure  5.7b  shows 
the  measurements  of  a  spray  formed  of  primary  droplets  and  doublets. 

The  theoretically  predicted  sizes  are  49  pm  and  62  pm.  The  sizes  actu¬ 
ally  measured  were  46  pm  and  57  pm  respectively.  Figure  5.7c  extends 
the  measurements  of  7b  to  the  presence  of  triplets.  The  theoretically 
predicted  sizes  in  this  case  are  49  pm,  62  pm  and  70  pm.  Note  that  the 
measured  diameters  of  the  doublets  and  triplets  are  related  to  the  pri- 
mary  droplets  by  2  '  and  3  respectively. 
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IMAX  DROPLET  SIZE  MEASUREMENTS 


FIGURE  5.7. 


Effect  of  beam  blockage  on  size  distribution 


The  effect  of  spray  blocking  the  laser  beams  before  they  cross  on 
the  size  distribution  measured  with  IMAX  was  explored.  A  size  range  of 
10  to  100  pm  was  used  In  this  case,  and  the  monodispersed  droplet  size 
was  73  pm. 

Figure  5.8a  shows  the  measurement  of  the  monodispersed  string  of 
droplets. 

Figure  5.8b  shows  similar  results  but  a  spray  is  blocking  the 
laser  beams.  Two  effects  can  be  noted:  the  peak  of  the  distribution 
dropped  to  69  pm  (5%)  and  the  spread  of  the  distribution  is  +5  pm,  -7  pm 
(+7%,  -10%),  to  the  1/e  .  It  should  be  noted  that  once  the  high  voltage 
is  corrected  to  account  for  the  beam  blockage,  the  broadening  should 
have  very  little  effect  in  the  distribution  of  a  polydispersed  spray. 

Spray  measurements  with  IMAX 

Measurements  of  a  spray  produced  by  a  pressure  nozzle  (Spraying 
Systems  TG0.3  at  50  psi)  were  conducted  with  IMAX.  They  were  made  at 
30  mm  from  the  nozzle  tip  and  at  two  radii:  0  and  10  mm.  The  results 
are  shown  on  Figures  5.9a  and  5.9b.  The  trend  shown  by  the  two  plots  is 
the  expected  one.  That  is,  there  are  more  small  droplets  in  the  middle 
of  the  spray  than  on  the  edge. 

In  order  to  test  the  resolution  of  the  system,  data  were  obtained 
using  three  different  size  ranges:  5  to  50  pm;  10  to  100  pm;  and  20  to 
200  pm.  This  is  one  of  the  most  difficult  self-consistency  tests 
imposed  on  any  technique,  and  most  available  techniques  will  show  a 
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55  100  10  55 

DIAMETER  (uni)  DIAMETER  (urn) 


NOZZLE:  SPRAYING  SYSTEMS  FULL  CONE  TG0.3 


shift  on  the  predicted  data.  IMAX  shows  excellent  matching  of  the  data 
in  the  overlapping  region  as  illustrated  in  Figures  5.9a  and  5.9b. 

The  corresponding  droplet  velocity  distributions  are  shown  on 
Figures  5.10a  and  5.10b.  It  can  be  observed  that  the  droplets  are  mov¬ 
ing  faster  at  the  edge  of  the  spray  than  in  the  center.  This  results 
from  the  large  momentum  and  terminal  velocity  associated  with  the  large 
droplets  present  at  the  edge  of  the  spray. 
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SPRAY  VELOCITY  DISTRIBUTION  WITH  IMAX 
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Figure  5.10b 
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